Extended Self-similarity in Kinetic Surface Roughening 
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We show from simulations that a limited mobility solid-on-solid model of kinetically rough surface 
growth exhibits extended self-similarity analogous to that found in fluid turbulence. The range over 
which scale-independent power-law behavior is observed is significantly enhanced if two correlation 
functions of different order, such as those representing two different moments of the difference in 
height between two points, are plotted against each other. This behavior, found in both one and two 
dimensions, suggests that the "relative" exponents may be more fundamental than the "absolute" 
ones. 
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Scale-invariant spatio-temporal behavior is observed 
in a wide variety of far-from-equilibrium systems. In 
analogy with the "universality" found in the equilibrium 
scaling behavior of systems near a second-order phase 
transition, it is interesting to enquire about the similari- 
ties between different non-equilibrium systems exhibiting 
scale-invariant behavior. In this paper, we point out a re- 
markable similarity between the scaling behavior of two 
well-known and extensively studied non-equilibrium sys- 
tems - turbulent fluids and growing interfaces. Krug |l]] 
discovered a similarity between the intermittent multi- 
scaling behavior of structure functions in strongly devel- 
oped turbulence ^|J|] and the scaling properties of cor- 
relation functions of height fluctuations in simple solid- 
on-solid models of kinetically rough epitaxial growth Q 
with limited surface mobility. The multiscaling proper- 
ties of these growth models have been subsequently inves- 
tigated in detail ]f| and a mechanism for this behavior 
has been proposed ||. In this paper, we demonstrate 
that the extended self-similarity (ESS) Jjj exhibited by 
the structure functions in fluid turbulence is also present 
in the behavior of correlation functions of height fluctua- 
tions in these growth models, and thereby establish that 
the analogy between deterministic turbulence in fluids 
and stochastically driven interface growth is remarkably 
deep. We emphasize that the ESS phenomenology in our 
discrete stochastic growth model is formally identical to 
that found in the intermittent fluid turbulence problem, 
establishing a precise one to one correspondence between 
these two seemingly completely different physical pro- 
cesses. While the exact reasons for this precise analogy 
between these two distinct problems remain unclear at 
this stage, we speculate that the existence of an infinite 
number of relevant (marginal) operators in both cases 
may be the underlying mathematical cause for this anal- 

ogy @. 

We begin by pointing out the analogy between fluid 



turbulence and surface growth. In fully developed turbu- 
lence, scaling behavior is observed in the inertial range 
r\ <C r <C L, where r is the length scale of interest, L is 
the outer integral scale at which energy is injected into 
the system and r\ is the inner dissipation scale. A mea- 
sure of the separation between the inner and outer scales 
is provided by the Reynolds number Re oc (-L/77) 4 / 3 . The 
g-th order longitudinal structure functions are defined as 

D,(r) = ([ti(x + r,t)-«(x > t)]«) 1 (1) 

where u(x, t) is the component of the velocity at the po- 
sition x at time t in the direction of the relative dis- 
placement r, and the brackets (. . .) represent a spatio- 
temporal average. These structure functions are believed 
to exhibit power-law scaling in the inertial range: 

D g (r) k r c ", i)<r<L. (2) 

The value of (3 is known to be exactly unity. The mea- 
sured values of Ca, q 7^ 3 differ appreciably from the Kol- 
mogorov result ||, Q = q/3, the deviation 5( q = ( q —q/3 
being positive for q < 3 and negative for q > 3. This is 
the phenomenon of multiscaling which is believed |^|,^| to 
arise from the strongly intermittent (violently fluctuat- 
ing) character of the local energy dissipation rate e(x, t) 
whose spatio-temporal fluctuations may be characterized 
by the quantities (e 9 (x, t))/(e(x, t)) q . These quantities 
are expected to exhibit a power-law dependence on the 
Reynolds number: 

(<*)/(*)* « [L/rjY" , (3) 

where fi q > for q > 1 . 

In models [[j],[l]^6) of growing interfaces, the role of the 
velocity field u(x, t) is played by the variable /i(x, t) that 
represents the height of a d-dimensional interface at point 
x at time t. In these models, the inner length scale is 
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the lattice spacing do which is usually set to unity. The 
outer scale is set by a correlation length £ that initially 
grows in time (£ oc t 1 ^ where z is the dynamic expo- 
nent) and eventually saturates (for t ^S> L z where L is 
the system size) to a value of the order of L. Scale- 
independent behavior of correlation functions is observed 
for length scales 1 <C r <C £. The correlation length £, 
therefore, plays the role of the Reynolds number. The 
height-difference correlation functions 



G q (r,t) = (\h(x + r,t) - h(x,t)\«) 1 /' 1 , 



(4) 



where (...) now represents a spatial average, are anal- 
ogous H to the structure functions defined in Eq.(0). 
These correlation functions exhibit the following scaling 
behavior: 



G q {r,t) 1 « r <£(*). 



(5) 



The observation |[|@| that C g is a decreasing function 
of g indicates that these correlation functions exhibit 
multiscaling similar to that found in turbulence. The 
quantity that is analogous to the energy dissipation rate 
e ~ (du/dx) 2 is the nearest-neighbor height difference 
s(x, t) = \h(x!,t) — /i(x, t)\ where x' is a nearest neighbor 
of the lattice point x. In analogy with Eq.(pl), different 
moments of s scale with different powers of the correla- 
tion length £: 



*M) = ([8(x,t)]*) i '*KZ(t) a < 



(6) 



where a q are found to increase with increasing q. 

In experimental and numerical studies of turbulence at 
relatively low Reynolds numbers, the range of values of 
r for which the power-law scaling of Eq. (j|) is observed is 
often very small and sometimes non-existent. ESS in tur- 
bulence refers to the fact that the size of the scaling re- 
gion is significantly enhanced if log \ D q (r)\, q ^ 3 is plot- 
ted against log |£>3(r) | . Since the result Z?3 (r) oc r in 
the inertial range is exact for the Navier-Stokes equation 
underlying fluid turbulence and log-log plots of |D 9 (r)| 
against |Z?3 (r) | show linear behavior over a substantially 
larger range than log-log plots of \D q (r)\ against r, ESS 
plots provide a convenient way of determining the val- 
ues of the multiscaling exponents f 9 . ESS has been used 
extensively Q during the last few years ]To|~|T2]| to an- 
alyze the data obtained from a variety of experiments 
and simulations. Substantial enhancement of the scaling 
region has also been observed in log-log plots of \G p (r)\ 
vs \G q (r)\ where p and q are any two unequal positive 
integers. The observation of ESS implies that the self- 
similarity of the velocity field extends beyond the con- 
ventionally defined inertial range. 

We have found that the correlation functions of height 
fluctuations in the Das Sarma-Tamborenea (DT) model 
Q of epitaxial growth in one and two dimensions ex- 
hibit properties that are very similar to the ESS described 
above. Krug's original work Q] on intermittent multiscal- 
ing in epitaxial growth was based on the one-dimensional 




FIG. 1. The height-difference correlation functions 
G q (r), q — 1 — 4 at time t = 10 as functions of the sepa- 
ration r for the Id DT model. The straight lines are best 
power-law fits of the data for r < 10. Inset: Rules of the 
DT growth model. Random deposition on the substrate is 
followed by nearest-neighbor diffusion if it is allowed by the 
local coordination, as indicated by the arrows. 



(Id) DT model. In this model, atoms are deposited ran- 
domly on a flat substrate under solid-on-solid condition. 
If a deposited atom has at least one lateral neighbor, 
it stays at that site. Otherwise, the atom moves to a 
nearest-neighbor lateral site if it can increase the number 
of lateral neighbors by doing so. If more than one such 
sites are available, then the atom moves to any one of 
them with equal probability. If no such nearest-neighbor 
site is available, the atom stays at the deposition site. 
These rules are illustrated in the inset of Fig.l. The DT 
model has been studied extensively jl],(|j|] by simulations. 

We first consider the scaling properties of the correla- 
tion functions G q (r) defined in Eq.(|j). Fig.l shows the 
results for G q (r), q = 1—4 at time t = 10 5 (in units 
of number of layers deposited) for the Id DT model with 
size L — 1000. The log-log plots of G q (r) against r clearly 
show linear behavior for small r, in agreement with the 
power-law form of Eq. (^]) . The values of the exponents 
Q' q obtained by fitting straight lines to the initial linear 
parts of the plots are: Ci = 0.88 ± 0.02, Q' 2 = 0.67 ± 0.02, 
C3 = 0.54 ± 0.02, and (' 4 = 0.45 ± 0.02. The size of the 
scaling region is, however, rather small (less than one 
decade in r). This is not surprising because the value of 
the correlation length £ oc t 1 ^, z ~ 4 is the order 

of 15 at t = 10°. The new and surprising result we find 
is that the data exhibit nearly perfect power-law scaling 
over the entire range 1 < r < 100 if G q (r), q ^ 1 is plot- 
ted against G\{r) (see Fig. 2). Since there is no "special" 
value of q in the growth model (in the sense that q = 3 is 
special in turbulence), we have chosen to plot the corre- 
lation functions for q ^ 1 against G\ - other choices lead 
to equally good scaling. It is clear from Fig. 2 that the 
range over which scaling behavior is observed is increased 
by more than an order of magnitude if one of the corre- 
lation functions is plotted against the other, rather than 
plotting them as functions of r. This model, therefore, 
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FIG. 2. Log-log plots of G q (r), q = 2 - 4 of Fig.l against 
Gi(r). The straight lines are best fits to power laws. Up- 
per inset: The local slope dlog[G4(r)]/dlog[Gi(r)] of the 
q — 4 plot in the main figure as a function of r (crosses), and 
same quantity for L = 100 in the saturation regime (squares) . 
Lower inset: The scaling functions g q (r), q= 1 — 4 (see text). 

exhibits ESS that is formally identical to that found in 
turbulence by Bcnzi et al . The extension of the scal- 
ing range is mostly towards large scales, which is similar 
to the behavior found in turbulence. 

The quality of the power-law scaling in the G q , q ^ 1 
versus G\ plots is illustrated in the upper inset of Fig. 2 
for q = 4, where we have shown the variation of the 
local derivative <f log(G4(r))/<f log(Gi(r)) with r. The 
observed variation is less than 10%, indicating that the 
power-law relation, G q (r,t) oc [Gi(r, t)] 0q , q ^ 1, pro- 
vides a very good description of the data over the full 
range of r. The values of the "relative" exponents 6 q , 
obtained from straight-line fits to the plots of Fig. 2, are: 
6*2 = 0.79 ± 0.02, 6» 3 = 0.64 ± 0.02, and 4 = 0.53 ± 0.02. 
These values are consistent with the expected result, 
6 q = Cg/Ci- The upper inset of Fig. 2 also shows the lo- 
cal slope of a log-log plot of G4 (r) versus Gi (r) obtained 
in the saturation state of samples with L = 100. The 
two sets of results are nearly indistiguishable, showing 
that the values of the exponents 9 q obtained from ESS 
plots are not sensitive to details such as sample size and 
the length of the simulation. Similar "universality" (i.e. 
insensitivity of the values of the relative exponents to de- 
tails such as the value of the Reynolds number and the 
flow geometry) has been observed in turbulence Jll]{i~2| • 

The observation of ESS implies a specific relation 
among the scaling functions [jlj that describe the behav- 
ior of the correlation functions G q (r,t). The dependence 
of G q (r,t) on r and t in the regime £ rj t x l z < L is 
expected to be described by the scaling form Jl]] 

G q {r,t)=t a «' z r^f q {r/t 1 l z ). (7) 

The occurrence of ESS is possible only if the scaling 
functions f q are related to one another by f^r/t 1 /*) — 
C q {t)\f{r/t 1 / z )]^i. This relation, when combined with 
Eq.(0), leads to the prediction that plots of the quan- 



FIG. 3. Log- log plots of G q (r), q — 2 — 4 against Gi(r) for 
the 2d DT model at time t = 10 3 . The straight lines are best 
fits to power laws. Inset: Log-log plots of G q (r), q — 1 — 4 as 
functions of the separation r for the same data. 

tity g q (r,t) = [G q (r,t)/G q (r — l,i)p^< against r should 
coincide for all values of q. The lower inset of Fig. 2 
shows plots of g q (r) (obtained from the data shown in 
Fig.l, using the values 0.870, 0.673, 0.540 and 0.452 for 
Q' q , q = 1 — 4, respectively) against r for q = 1 — 4. The 
data for different q collapse nicely to the same curve, 
confirming the occurrence of ESS in this system. A sim- 
ilar description of ESS in turbulence in terms of a q- 
independent scaling function is provided in Ref. |i"o| . 

The behavior described above is also found in the 2d 
DT model, as shown in Fig. 3, where we have plotted 
G q {r), q = 1 — 4 against r (inset), and G q {r), q — 2 — 4 
against G\(r) (main part) on log-log scales. The data 
shown correspond to |5j t = 10 3 for 500x500 samples. We 
again find an extension of the scaling region by more than 
one decade. The values of the exponents obtained from 
fits to the data are: 9 2 = 0.77 ± 0.02 , 9 3 = 0.60 ± 0.02, 
6> 4 = 0.48 ± 0.02. As before, the relation 6 q = C^/Ci is 
satisfied within the error bars. 

We have also found ESS in the dependence of the quan- 
tities o q denned in Eq. (|^) on time t (in the growth regime) 
or the sample size L (in the saturation regime) . The inset 
of Fig. 4 shows our data for the dependence of a q , q = 1—4 
at saturation on L in the 2d DT model. The log-log 
plots clearly show a downward curvature (similar cur- 
vature is observed in log-log plots of a q versus t in 
the growth regime) that becomes more pronounced as L 
is increased. In contrast, the ESS plots of a q , q — 2 — 4 
against cr±, shown in Fig. 4, do not exhibit any such curva- 
ture. The values of the "relative" exponents ip q , defined 
by o~ q oc (ai)^ q , obtained from straight-line fits to the 
the ESS plots are: ip 2 ^ 1.2, -03 - 1-5, V>4 ^ 1.9- 

The results shown in Fig. 4 and qualitatively similar 
behavior found in the Id DT model [l]|| clearly indicate 
that the DT model does not exhibit true asymptotic mul- 
tiscaling. This conclusion is corroborated by recent stud- 
ies H which suggest that the approximate multiscaling 
observed in these models is a non-universal and extremely 
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FIG. 4. Log- log plots of a q (L), q — 2 — 4 against <Ji(L) 
for the 2d DT model at saturation. The straight lines are 
best fits to power laws. Inset: Log-log plots of a q , q — 1 — 4 
as functions of the system size L for the 2d DT model at 
saturation. 

slow transient or crossover arising from a non-linear in- 
stability in the discretized version of the underlying con- 
tinuum growth equation. The ESS found in this paper 
shows that the slow crossover responsible for the approx- 
imate multiscaling behavior affects correlation functions 
of different order in exactly the same way. In this picture, 
the multiscaling exponents (,' q and a q are, at best, effec- 
tive ones describing the behavior over a limited range of 
length and time scales. The occurrence of ESS implies 
that "relative" exponents such as 9 q and ij) q are, in some 
sense, more fundamental than the "absolute" ones, Cf q 
and a q . While it would be premature to suggest that a 
similar transient description applies also to multiscaling 
in turbulence, we note that the possibility that the in- 
termittency corrections 6£ q are finite Reynolds number 
effects which would vanish in the Re — > oo asymptotic 
limit has received considerable attention in the re- 
cent literature. Also, there is some numerical evidence 
|L5j indicating that the "relative" exponents are more 
universal than the "absolute" ones in turbulence. 

Our finding that ESS may occur in problems (i.e. in 
the kinetic surface roughening of the DT model) very dif- 
ferent from the fully developed homogeneous turbulence 
problem (where ESS was originally discovered 0) is po- 
tentially significant, and may eventually provide a clue 
to its understanding. Currently the ESS phenomenon 
(both in surface growth and turbulence) remains an in- 
teresting empirical fact without any rigorous theoretical 
understanding. It has recently been shown || that mul- 
tiscaling in the DT model shares a substantial common 
phenomenology with that in fully developed turbulence, 
and the intermittent multiscaling behavior in these very 
different problems may arise || from the existence of an 
infinite number of marginal operators and an associated 
near-singularity in both problems. We speculate that 
the eventual theoretical understanding of ESS will de- 
pend on a more detailed understanding of the roles that 
the infinite number of relevant operators and the near- 



singularities play in homogeneous turbulence and kinetic 
surface roughening. It would obviously be very inter- 
esting to explore further the analogy between kinetically 
rough interfaces and fluid turbulence. 
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